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Abstract 

A physics-based methodology is described to predict jet- 
mixing noise due to small-scale turbulence. Both self- and 
shear-noise source terms of Lilley’s equation are modeled 
and the far-field aerodynamic noise is expressed as an 
integral over the jet volume of the source multiplied by an 
appropriate Green’s function which accounts for source 
convection and mean-flow refraction. Our primary interest 
here is to include transverse gradients of the mean density in 
the source modeling. It is shown that, in addition to the 
usual quadrupole type sources which scale to the fourth- 
power of the acoustic wave number, additional dipole and 
monopole sources are present that scale to lower powers of 
wave number. Various two-point correlations are modeled 
and an approximate solution to noise spectra due to multi- 
pole sources of various orders is developed. Mean flow and 
turbulence information is provided through RANS-ks 
solution. Numerical results are presented for a subsonic jet at 
a range of temperatures and Mach numbers. Predictions 
indicated a decrease in high frequency noise with added 
heat, while changes in the low frequency noise depend on jet 
velocity and observer angle. 

1. Introduction 

Development of a reliable and robust jet noise 
prediction method is instrumental in the design of jet noise 
abatement tools. Short of Computational Aeroacoustics 
(CAA) which is arguably the hope of the future [1, 2], most 
near-term approaches to the jet noise problem start with 
some derivative of Lighthill’s theory in order to describe the 
sources of aerodynamic noise as well as the propagation 
effects (e.g. chapter 6, Ref 3). Pilon and Morris [4] 
recently developed a semi-analytical method for jet noise 
prediction that employs scaling laws based on acoustic 
analogy to characterize the quadrupole and dipole sources 
and uses a simplified ray-acoustic approximation to 
propagate the sound to the far field. Tam et aL [5] 
concluded that a two-component spectra is needed to 
describe the noise generation due to both small-scale 
turbulence and instability waves in the jet flow. They 
maintain that instability waves generate sound at subsonic as 
well as supersonic speeds, although they become 
increasingly more efficient at high supersonic speeds. 
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More recently, Tam and Auriault [6] developed a 
semi-empirical theory to predict the small-scale turbulence 
noise based on the kinetic theory of gases. They argue that 
the source of small-scale turbulence noise is the kinetic 
energy of fine-scale turbulence and use a Reynolds 
Averaged Navier-Stokes (RANS) solution with a ke 
turbulence model to obtain the source strength as well as 
the time- and length scales of small eddies. 

Although, in principle, direct numerical 
simulation (DNS) may be used to describe the acoustic 
near field [7,8], all scales of motion cannot be resolved for 
the high Reynolds-number jets of practical interest due to 
increasing demand on grid resolution. It is argued that at 
supersonic speeds, large-scale structures dominate the 
initial noise-producing region of the jet. Consequently, 
large eddy simulation (LES) has been advocated [9,10] as 
a useful and practical tool to predict the acoustically 
efficient larger scales while modeling the smaller subgrid- 
scale turbulence. 

In the present study, we employ a physics-based 
modeling ^proach to describe the soimd sources of small- 
scale tuibulence. Our primary concerns are the so-called 
self- and shear-noise sources of Lilley’s equations, 

commonly represented as pDd^{u.Uj)l dx^dxj and 

-2p{dU ! dr)d^{u^u^!dxpx^ respectively. A great deal of 
simplification in modeling of a two-point correlation may 
be achieved by assuming that the mean density p is 
constant (i.e, quasi-incompressible turbulence), and by 
considering turbulence as isotropic. Reference [11] details 
the derivation of various quadrupole type sources that 
emerge from permutations of a two-point velocity 
correlation using both isotropic and axisymmetric 
turbulence models. Here, we assume that die mean flow is 
locally parallel and its density is a function of the spanwise 
coordinate r. The radial gradient of the mean density (and 
sound speed) produces additional sources within both the 
self- and shear-noise terms. Approximate expressions for 
noise spectra from various source types is given and the 
significance of newly derived sources and their 
contribution to noise of heated jets is examined through a 
numerical simulation of a simple roimd jet. 
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On the experimental side, the effect of jet density 
on aerodynamic mixing noise from subsonic and fiilly 
expanded supersonic jets has been studied by a number of 


investigators. Hoch et al [12] studied the effect of density 
on noise by varying the jet temperature at a range of 
velocities from 1 50 to 800 m/s. Their far-field 
measxirements show the expected reduction in noise with 
reducing density at high velocities and the existence of a low 
velocity regime in which the reduction of density by heating 
the jet increases the acoustic power. A more detailed 
experimental study of the effect of jet density on noise 
spectra is given by Tanna [13]. According to his findings, at 
high Strouhal numbers, mixing noise decreases with heating 
the jet throughout the velocity range. At low source Strouhal 
numbers, the intensity increases with heating at subsonic 
acoustic Mach numbers (i//a^) while it decreases with 
heating at supersonic Mach numbers. We intend to 
compare our predictions with these observations using jets 
numerically calculated at inflow boimdary conditions that 
could closely simulate the experiment. 

For the most part, our predictions agree with data 
and indicate a reduction in high frequency noise with 
increasing jet temperature. At intermediate and low subsonic 
speeds, heating enhances the low frequency noise at smaller 
angle from jet inlet. In addition, predictions indicate that 
density gradients of the mean flow contribute primarily in 
the context of the dipole sources of self noise and radiate 
mainly at 90®, 

The paper begins with the formulation of the 
governing equations. In sec 2.1 we briefly describe Balsams 
Green’s fimction. Modeling of two-point correlations and 
sound spectral intensity due to various terms of self and 
shear noise are described in section 2.2. Sections 3 and 4 
describe implementation of the prediction methodology and 
the final summary. 

2. Spectral Density of Sound as Related to Source 
and Green’s Function 

We assume float die mean flow is locally parallel and float 
local daisity and sound speed are functions of the radial coordinate 
ronly 

U = U{r), p = p{r), a = a{r). (1) 


whoe p is the acoustic j^essure, daiote the cylindrical 

coordinates, D-dldt^Udldx^ is the amvective derivative, 
V^is the Laplacian c^oerator in a cylindrical cocHxiinate system 
and W, is the fluctuating part of the velocity. The terns on the 
right-hand side of Eq. (2) are usually referred to as self- and 
shear-noise source terns, respectively. 

2.1 Moving Reference Frame (Balsa’s Solution) 

To axount for both source convecticm and the 
refraction of sound in the jxesence of the mean flow, a closed 
form sohiticm to Lilley’s equation in flie hi^ frequency limit was 
discussed by Balsa in References 14, 15 and 16. Balsa’s sohiticHi 
is exfxessed in a coordinate system T moving wifli flie 
tmbulait eddies, 

x^=x-Uj, V = U-U,, D^dldt^Vdldx:,. (3) 



Fig. 1 Moving frame 

The ccMivectioo velocity is in thexj direction and is usually 

some fraction of jet exit velocity. Here, Q refers to the source 
frequency in a moving reference frame, Le., a frame moving with 
a ccmvecting eddy. A two-point space-4ime correlation is 

denoted as R(y,^, r) where ^ and r represent the ^>atial and 

ternporalsq3arati(Xibetwe«i points j?j=y-^and 

Experimentally, the coiTelation R(y^ r) in a jet flow describes 
a fluctuating pattern in a mo viiig frame and is expressed as 


Lilley’s equatm for an inviscid flow linearized about 
imidiectional transversely sheared mean flow is 


Lip; U,x,)^\ D^p -DV^p--^ (log )D 
a dr 


dr 


a 


dr 

dxfdxj ^ dr dxpx^ 

SeIf-Nois€ Shear-Noise 


( 2 ) 


= l-Uj. 

Far-field spectral density beccnnes 

7(jc',Q) = JjG’(r,j;-|/2.Q)G(x',J5 4|/2,Q) 

yi. 

(4) 

wliere * indicates a omiplex cMijugate and iy, , ^2) is 
the ^)ectra] density of a two-point space-time ccnrelation in a 
moving fiame 
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—00 

A Doppler fector relates source frequency Q to the observer 
frequency cu. 

For teevity, we shall suppress subscr^ m in writing 

source correlatictti functiais and use ^ as a sg)arati<xi vector with 
respect to a moving frame in the remaining discussions. It turns 
out fliat the Greoi’s fimcticHi solution to Eq. (3) depends on the 
behavior of a so-called shielding functicm g in frie vicinity of flie 
source. The shielding function is defined as 


g\r,0) 


(1- Af cos^)"(q, ! of -cos^ e 
(1- M^cos^)^ 


( 6 ) 


Ibe above Green’s function was obtained for a &r-field observer 
(Le. R! D»\). When the shielding function crosses zero 
(xify race, Le. g(r^,ff) = 0, and if die source is located at 
r^<r^, die shielding fimcdon ^(r,^ becomes negative fw 
r<r^. The Green’s fimctKHi is dieiejqxessedas'^^ 


G(x , j/,Q) = 


i aola^ 

AmJ(R (1 - Mg cosfff 




y^ofoj 


xe e (7c) 

whaie 

f=-g\ C = \/ir,e)dr. (7d) 


where M=^(r)/a^, ^ source cmwecdoa Madi 

number, 0 is tiie polar angle fixm downstream jet axis and 
subsapt 00 refers to ambiait The zeros of g^, denoted as ,are 

called the turning points. Shielding of radiated sound d^)oids on 
flie number of turning points in the vicinity of a source. When 
there are no turning points (i.e, g^ >0); the Green’s functicxi 
solution to Eq. (3) ( referred to as an unshielded solution) is 


G(x\y,Q) = 


i Op la. C ^ 


xe 


g-g» )] 




(la) 



Fig. 3 One turning point problem. 



Fig. 2 No turning point iMoblan. 
and 

^=\[g(r,ff)dr, A: = Q/a„, 

g^ = g{r -^(0,0) = sin0/(l - Mg cos0) . 


Thus /^is positive wiien r <r^j (Fig. 3) and frxxn (7d) it is seen 
that the last exponential factcH* results in noise attmuation ^\ilich 
becomes stronger the deeper the source is embedded in the 
negative region ofg^fh^ and the higher its fiequaicy. The case 
with two turning points is also discussed in flie same references. 
In general, flie shielding of noise depends on flie proximity of the 
source with respect to the turning point r^as well as the number 
of turning points. 

22 Noise Sources 

The prediction of Jet noise due to small-scale 
turbulence may be carried out according to Eq. (4) once flie 
spectral daisity of a two-point correlaticHi associated with various 
source terms of Lilley’s equation are known. Any realistic 
c(Hnputati(xi of Eq. (4) requires SOTie sort of acceptable modeling 
to Militate a closed-fam int^ration with respect to the 

separation vector If the time history of source fluctuatiOTS is 

( 7 b) not known, additional modeling related to the temporal fector of 
a two-point correlaticn beoMnes necessary in order to e\^uate 
the source spectral density. In fliis secticm, we discuss flie two 
relatively mcxe significant sources of Lilley’s equation, namely 
self- and shear-noise terms as given in Eq. (2). 
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2^.1 Self Noise 


To assess die fer-field mean-square pressure due to the 
self-noise tenn onfy, we write (2) in a moving coordinate 


L{p„^-,v,x,) = r\ 


dx\dx'j J 


( 8 ) 


where density has moved to die right of operator 
Z) = 9/5/ -I- F5/cbc’ assuming fliat flow is locally parallel and 
that density fluctuaticws are small so is the mean density at die 

source. We now write Lfll^’s Green’s function S for a source of 

typeZ){e"'^9(x -x')} 


L{Se'^ ;V,x\) =D{e'‘^S{x)S{q>- q>, )5ij- yr] . (9) 
It is easily shown that S is related to G 


In (13b), the integrand of tiie source term 7^.^ is the familiar 
fourth-order g>ace-time correlation. The velocity variables with 
and without prime are evaluated at two points andy^ separated 

by vectOT ^ and time-delay r . The factor multiplying the 
source intensity in (13a) defines a directivity factor that accounts 
for effects such as source convection and mean-flow refraction. 

The volume integration in (13a) needs to include the 
most energetic parts of fiie jet When the flow is axisymmetric, 
the directivity factor may be avoaged azimuthally with respect to 
source and observer circumferential angles to obtain a ring- 
source directivity fector Subsequently, jet volume 

integration will be limited to radial and axial coordinates 


y\ ^ 


(14a) 

(14b) 


1-M^cos^ 

The acoustic pressure due to die above source and Green’s function 
becomes 

d\u,Uj) 


= J J S{x ,f,y,t,)p^ ^ 


'i y 


( 11 ) 


where S{x inverse Fourier transform of 

S{x , y, Q)e'^ . Integrating by pats widi suitable restrictions 
on S and its derivatives at large values of y^ one finds 


Appendix A gives ^jpx)ximate expressions for S ,j 
upon neglecting /J^_^__and M , . Directivity factors are 

now calculated by placing afp-c^>riate variables S .j into Eq. 
(14a). Of die 81 canpements, it turns out dial some are equal and 
some are identically zero. Whai the shielding function g^(i') is 
positive for values of radial cocxdinate r fixim source point to 
the djserver point, die solution is refered to as unshielded and 
die corresponding directivity factors relevant to our calculations 
becOTie 




y 




(12) 


The compactness conditions requires Q/ / to be small 
(/ is a typical eddy laigth scale). Since fiequency O is of die order 
of v//(v is die R.M.S. of the velocity fluctuations) then the 
compactness ccBiditkMi should be true if v / < 2 ^ is small. 
Consequendy, in fine-grain turbulence, the eddy lengdi scale is 
expected to be much shorter dian the wavelengh of fte acoustic 
disturbances (/« X). If diis is true, then mean square acoustic 
disturbances in die fer field may be expressed as a superposition of 
disturbances from independently correlated volume elements in a 

jet 

Pl,/ix,0.) = (x,j,Q)S*,^ I.j„(y,ny^ , (13a) 

y 

where 

-KO 

lj„(y,Cl) = J J {u,Ujyu,u,y^^dTd^ . (13b) 


1,,ia|2 COS0 


«2222 =t|‘^ 


_ 1 ^2 
^2233 ~ 4 


- Me' 

<^2323 ” T F* 


3 

—(erf +Oj +al)-2ap^-a^a^+2ap^ + 

b^,+bl-b,b^ 

^ (af -t- Oj -I- a' ) - 2a,flj -i- 2a^a.^ - - 

^(a^+al +a^) + a^a,+hf +b\ +b,b^ 

(15a) 
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Factors and have a dimensicm of and depend on flie 
mean flowfTqjaties as 

«2=N, 03 =— (T,-^), 

2g„ 

b,=-2kgj,, (15b) 

Variables aj, ci^,b^ 2 nd become identically zero wiiengradiait 
terms/? „ and M areneglected(see Appendix A), and hence 
the directivity factcMS (15a) reduce to tiiose given in [11]. 

When g^{r) becomes negative at the source, a 
shielding zeme exists. The amount of shielding depends on the 
proximity of source wMi respect to the turning point as well as 
die number of turning points. This dependence is accounted for by 
multiplying directivity fectors (15a) by a shielding coefl5cient . 

For example v^en there is only one turning point , the shielding 
coeifficient becomes 

(15c) 

Eq. (15a) is now written for a unit volxime ring widiin a turbulent 
jet 


Fuidier, using the usual assun^cxis we relate die axial length 
scale Zj to the turbulence kinetic aaagy and its diss^iation rate c 

as Z, /s. hi addMcffli, the source chaiacteistic 

fiequCTicy (which is die inverse of charactaistic time-del^ 

) is written as - s/k . Time-delay fimcticxi of a two- 
point correlation was selected as 

g(r) = exp(--^(a/2)^+ cr = 0.8. Making 

die necessary substitutions, we findfron (17) 


• 13/2 


(18) 

where ZT] is a modified Bessel functicxiofthe second kind and of 
ordereme. A I>^lerfktCH* relates die source fiiequencyQ to 
the obsaver fiequency o as 

Q = 0)^(1 - cosfff + » 

where is an empirical ccxistant selected as 0.50. 

2.23 Shear Noise 

The sound field due to shear^ioise terns of Lilley’s Eq. 
(seeEq.2) is given using the Greoi’s fimedon of Eq. (7) 






^\ + {arJTf 3 


= +4a„j3/„j2 + 8^1212^212 

^^ 2222^2222 *^^^2233*^2233 “^^^2323*^2323. 


233 Source Correlations Tensors 


(16) 


Self noise is evaluated by integrating (16) for die most energetic 
parts of a turbulent jet Source cortelaticHi tensOTS /^.^(5/,Q) are 

usually calculated using a |hysics-based modeling ^^pxiach. As 
was done in Ref [1 1], we assume that turbulence is axisymmetric 
about die direction of the mean flow. Various source cOTrelaricn 
tensors are simplified and written in terms of the axial compcxi^t 
(see Appendix B) 


= (17) 

Here g(r) is the temporal part of a two-point corcelation u^Uj , 
Zj and Z 2 are axial and transverse length scales of an 
axisymmetric turbulence and A = Z 2 / Zj . Parameter 

j3=(\-ul/u^)is assumed known a priori, and the axial 
COTiponent of turbulence inteisity is expressed as 
— 2 2 

= -jX'/(l-— ^) where vis the turbulence kinetic energy. 


Pshear(x,t) = ~2j J G(x )p 

h y 


dUdi^ju^Uj) 
dr dy^dyj 


dt,dy 


(19) 


It is implied that subscr^t 2 on the source tern retes to a 
transverse direction, whidi when expanded is rqilaced with 


+ Integration Ity parts results in 




h y 


In writing (20), we have used the parallel flow assumption 

dpldy ^-^ , and dWdy^ 0 . Let G denote the quantity inside 
the bracket in the last equation. Upcxi recognizing that 
dGIdy^ =/^cos^/(l- A/^cos^), we write 


^ dUdG 

Cr = , 

“ dr 

+^r+JT 


( 21 ) 


pLar(x,Oi) = J J b,jl^,^j{y,^)drdy, 


i'l ' 
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Upon neglecting the second derivative U ^ , expressions for G ^ 
become 




^ cos^ ' 




\- M^cosO 


^ 0 ^ 


= ik 

-g„cos<p 

G- 

Tj cos^p„ 



-g„smq> 

J 




and 


p , M, COS0 

Y _ I n >^o 

^ 2/?^ 1— Af^cos^ 



(22a) 


(22b) 


Using (22a) in (21), shear noise due to a unit volume ring source 
becOTfies 




Fig. 3 RC nozzle ge(Mnetiy. 

A tc^ of four run points are diown here. Tabic la 
shows that by lowering the nozzle pressure ratio (NPR) fiom 1.4 
to 1.27, and at die same time increasing iqjstream total 
temperature TJT^ fixm 1.55 to 227, the jet exit velocity may be 
kept ccMistant at Gy/a^=0.83 while exit temperature 7) is 
increased fixim 1.42 to 2.13. Similar run points are shown in 
Table lb for a low subsonic jet velocity of V =0.52 


P shear (^> ^5 3 ^) ^(^1 1 ^2\2\ ^ 12 ^2232 ^ 13^^2323 ) ’ (^^ 


and the corresponding unshielded directivity fectos are 




COS0 


\-M^cose dr 




COS0 


1-ATcos^ dr 


^13 ” ^12 * 


(24) 


For the shielcfed case, (Le.when g^{r^O 

shielding factor fd^j as discussed earli©- should mult^ly the above 

directivity fectors. 

Whengradients p and Af are neglected, it is 

readily seen that die directivity fecttss of shear noise also reduce to 
those given in Ref. [1 1]. 


3. Numerical Results 


Table la Inflow conditioiis for exit velocity of r^oo=0.83 



Inflow Conditions 
(inpiut) 

Exit Ccffiditions 
(calculated) 

Case 

NPR 

TJT^ 


TJT^ 

J 

1 

1.40 

155 

0.83 

1.42 

2 

127 

221 

0.83 

2.13 


Table lb Inflow conditioiis for exit velocity of U/a^=QSl 



Inflow ConditicHis 
(inp>ut) 

ExitCondMOTS 

(calculated) 

mu 


TIT 

Ujla^ 

TJT^ 



155 

052 

1.50 

4 

1.095 

221 


221 


Predicted radial frofiles fcx* mean velocity, ten^ierature, 
and density are shown (Fig. 4) at die intennediate jet exit velocity 
of Ujla^ =0.83 and at four streamwise locadcms akmg the jet 

axis; I^ce the effect of heating on jet noise is examined by 
eliminating velocity dq^endence frOTi calculations. 


A single stream Round-O)nvergent nozde (RQ widi an 
exit diamete of Z>=12,13 cm. is selected fcr numerical simuIatMXi 
(Fig. 3). RANS calculations are performed with the NPARC 
Navier-Stokes solver [17] and Qdm's k-e turbulence model [18]. 
Inflow ccMiditions arc specified at inlet (4.33D fiom nozzle exit), 
and downstream conditicxis are provided at 20D fiom exit plane. 
In a spanwise directicMi, the grid extraids 8 . 1 1 D frcMH the centerline, 
and a fiee-type bwmdaiy conditicm is specified Our primary 
objective hoe is to examine die effect of temperature, or 
equivalaitly the mean density gradient on radiated jet noise. A set 
of run points for die CFD predictions are carefiiUy selected with die 
intenticm of changing die jet toiperature while maintaining its 
velocity. 


Noise predictions were pierfcxmed at subsonic 
ccmditions of Tables la and lb, and oxnpiared with expierimental 
data [13], Measurements of Tanna [13] wctb conducted for a 
convergent nozzle with exit diameter of 5.08 on and cm an arc at 
72 diametes fixHn exit, AWiough the Stalls ofthe geometry and 
run px>ints for our numerical simulaticm are not exactly die same 
as die expieriment, for all practical purpxises a similar piarametric 
stu^ is conducted and a qualitative a»np>arison is justifiable. 
For example, measurements indicate diat at hi^ fiequency, 
mixing noise decreases with heating the jet thrcHighout the jet 
velocity (see Fig. 5). Low fiequency noise increases with heating 
at subsonic acoustic Mach numbers Uj/a^(F\g. 5), while it 

decreases widi heating at stposonic Mach numbers (not shown 
here). Measured directivity ^ig. 6) indicates an ovoaIlreducti<Hi 
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in noise with headng at high and intamediate velocity and an 
increase in noise at low subsonic Mach number of U Ja =0.50. 

j « 

Figures 7 and 8 show flie effect of heating on the 
predicted spectra at intermediate and low jet velocities of 
Uj/a^ =0.83 and 0.52. The spectra are ccHnputed at /M>=125. 

Clearly the temperature gap between the jets in our simulaticHi 
(Table l)isnotaswideasinthee?q5erim«it,yetanincreaseinlow 
fiequency noise and a reduction in high frequency noise is still 
evident The integrated spectra (Fig. 9) shows a reduction in SPL 
wifli increasing temp^ature at all angles. At low Mach number of 
Fig,10, the directivity facte does not exactly follow the trend 
suggested by data. Comparingfigures7and8 withdataofFig. 5, 
it is observed that the increase in low frequency noise in our 
simulation in not as large as suggested by data. However, the 
density gradioit tenns jHOvide improvements at the low end of the 
spectra and suggest a trend more compatible wMi measuranait 
Figures 1 1 and 12 demonstrate die effect of heating on noise 
directivity at specified Strouhal numbers (St==fi)/Uj). They show 
that at low frequency, heating results in a relatively small increase 
in noise intensity at mid angles. At high fiequency, addition of heat 
reduces the noise and increasingly more so at a higher jet velocity. 
Fig. 13 examines the contributions from muM-pole axnpcments of 
self and shear noise. Noise components Q, D and M denote the 
quadrupole, dipole and raonqx)Ie predictions, which are 
proportional to k\ and respectively. When density gradient 
terms are neglected, only die quadn^le contributions are pesent 
Outside that, the new sources ^ipear primarily at the low end of 
die spectra. 


Appendix A 

Various derivatives of 5 are needed with respect to 
source coordinates . We use the equation of state fcH* an ideal 


gas and exp-ess the sound speed as a\-Pylp ^ , where P is a 

constant static pressure; hence S is poportional to . In 
addition, die radial gradients of die Mach number and sound 
speed appearing in die Green’s function (see Eq, 7) are accounted 
for to the first derivative. A transformation is made to a polar- 
cylindrical cOOTidinate system 
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Further increase in jet speed in the present numerical 
simulation requires a shodk-fiee convergent-dvergent nozzles 
(CD). An increase in die jet temperature at a given supersonic 
speed requires anew design Mach number and hence a new nozzle 
geometry. Measurements [13] indicate that noise level should drop 
with heat at die entire range of die fiequency ^pectnim. 


4. Conclusions 

This stucty was motivated by a desire to examine the 
effect ofheating (Ml jet noise spectra In additi<Mi to the mean flow 
refiaction effects that directly influence the directivity of jet noise, 
new sources related to the radial gradients of mean density and 
sound speed were introduced A fhysics-based modeling 
approach was develc^ped to proximate self- and shear-noise 
sources of Lilley’s equations. Turbulent vel(x:ity fluctuations were 
assumed axisymmetric about the direction of mean flow. The 
present model enhanced die prediction capabilities of jet mixing 
noise due to small-scale turbulence. It was suggested outside 
die usual quadrupoles, c<Mitributions firm newly calculated sources 
were mainly due to dipole sources of self noise and at the low end 
of the fi^equency spectrum. 
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and N is given as 
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Once die ratio a ^ is replaced with -p ^ /(2pJ and radial 
gradients of variables Tj and T 2 become 
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Appendix B 


Various source correlation tensor are related to 


the axial tensor I,, 


^2222 

*" ^3333 “ 

Cl 

^U22 

” A 133 ” 

C I 
'-'2^1111, 

^2233 

-Cl 

^3‘*1UI, 


^313 

” ^1212 ~ 

C I 

A 323 

= C I 

. 


m 


For an axisymmetric turbulence coefficients Q are given as 
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Fig. 4 Comparison of mean flow profiles for conditions of Table la and at indicated 
streamwise locations. 
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Fig. 13a Predicted multi-pole ccHiqxments of self no^ 
at q^erating conditions of Table.la case 2. 


Fig. 13b- Predictedmulti-polecOTipcxients of diear noise 
at q)CTating conditions of Table. 1 a case 2. 
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